Dimensionality effects in restricted bosonic and fermionic systems 
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The phenomenon of Bose-like condensation, the continu- 
ous change of the dimensionality of the particle distribution 
as a consequence of freezing out of one or more degrees of free- 
dom in the limit of low particle density, is investigated theo- 
retically in the case of closed systems of massive bosons and 
fermions, described by general single-particle hamiltonians. 
This phenomenon is similar for both types of particles and, 
for some energy spectra, exhibits features specific to multiple- 
step Bose-Einstein condensation, for instance the appearance 
of maxima in the specific heat. In the case of fermions, as the 
particle density increases, another phenomenon is also ob- 
served. For certain types of single particle hamiltonians, the 
specific heat is approaching asymptotically a divergent behav- 
ior at zero temperature, as the Fermi energy ef is converging 
towards any value from an infinite discrete set of energies: 
{ei}i>i- If £f = e», for any i, the specific heat is divergent at 
T = just in infinite systems, whereas for any finite system 
the specific heat approaches zero at low enough temperatures. 
The results are particularized for particles trapped inside par- 
allelepipedic boxes and harmonic potentials. 

PACS numbers: 05.30. Ch, 64.90. +b, 05.30.Fk, 05.30.Jp 



I. INTRODUCTION 

Because of the advances in nanotechnology it has be- 
come possible to use very small structures in a broad 
range of applications. The importance of these applica- 
tions and the fact that the physical properties of such 
structures could be very different from those of bulk ma- 
terials, make the theoretical and experimental investiga- 
tions very useful in this area. 

The experimental findings in Q motivated us to calcu- 
late the thermal properties of ultrathin dielectric mem- 
branes or wires by splitting the phonon spectra into 
discrete and a continuous parts This framework 

implies crossovers between different phonon gas distri- 
butions, reflected, for example, in the exponent of the 
temperature dependence of the specific heat or heat con- 
ductivity. For example in a membrane, as the tempera- 
ture drops, the population of the phonon modes parallel 
to the surfaces [which we shall call the two-dimensional 
ground state (2D gs)] becomes dominant, and the three- 
dimensional (3D) phonon gas distribution changes into a 
two-dimensional (2D) one 0. The macroscopic popula- 
tion of the 2D gs [or one-dimensional ground state (ID gs) 
in the case of a wire ||] and the qualitative differences 



between phonon gas distributions with various dimen- 
sions enabled us to make the analogy with the multiple- 
step Bose-Einstein condensation (BEC) |||| and to call 
this phenomenon Bose-like condensation (BLC). Yet, the 
number of phonons changes with temperature and fea- 
tures like maxima of the specific heat (cy) observed in 
the case of BEC can not be seen in the case of a phonon 
gas undergoing BLC. 

The first purpose of this paper is to extend the pre- 
vious work reported in Refs. [g|| and to describe BLC 
in systems of massive bosons and fermions. This will 
be done in Section O. The mathematical technique used 
here is a straightforward extension of the one introduced 
by Pathria and Greenspoon in Ref . H . Nevertheless, the 
analytical approximations used there are not appropriate 
for our case. Therefore, after obtaining general expres- 
sions, we make numerical calculations to give concrete 
examples of BLC and to observe the behavior of the spe- 
cific heat during the transition. The phenomenon occurs 
at low particle densities (this will be made more clear in 
Section |l]) and is specific to both bosons and fermions. 
At low temperatures, the number of massive particles in 
a closed system can be considered to be constant. The 
conservation of the particle number will allow us to ob- 
serve resemblances with the BEC, like, in some cases, 
maxima of the specific heat (cy max ) a t the condensation 
temperature. Anyway, the signature of BLC, as seen in 
the temperature dependence of cy, is more complex and 
depends on the energy spectrum. 

A consequence of the third law of thermodynamics is 
that the specific heat of any thermodynamical system 
should vanish at zero temperature. Li et al. showed in 
Ref. Q that the heat capacity of a Fermi gas, confined in 
an external potential of quite general form, and for any 
space dimension, has the asymptotic behavior cy oc T 
at low temperatures (where T is the temperature of the 
system). This is for the case of a continuous energy spec- 
trum. In contrast to this we show in Section III that the 



specific heat of a Fermi gas with a single-particle hamil- 
tonian of the form H = H c + H^, with H c having a 
(quasi)continuous spectrum e c S [0, oo) and having 
the discrete eigenvalues e,, i = 0,1,..., may approach, 
depending on the density of the energy levels of H c , di- 
vergent behavior at temperature T = K as the Fermi 
energy converges to ej, for any i > 1. In such a case, if 
the spectrum of H c is continuous, then the specific heat 
diverges at T — and ep = e,-, for any i > 1. However, 
in any finite system the energy spectrum is discrete, so 
the specific heat approaches zero if we go at low enough 
temperatures and the third law of thermodynamics is not 
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violated. 

Ultrathin (semi)conducting membranes and wires, 
nowadays widely used in mesoscopic applications, atoms 
in very anisotropic harmonic traps, wires or constrictions 
defined in 2D electronic gasses are just a few examples 
of systems where the phenomena presented here could be 
observed. Also, they could provide an understanding of 
the behavior of very thin liquid He films. 



II. BOSE-LIKE CONDENSATION 

The BEC in cuboidal boxes with small dimensions 
drew a lot of attention many years ago, in the beginning 
in connection with very thin films of liquid He ||§,§-jLlj . 
It is now well known that, as the dimensions of the box 
are reduced, at constant density, the cusp-like maximum 
of cy is rounded off and the condensation temperature (in 
this situation taken as the temperature corresponding to 
the maximum) increases with respect of the bulk value. 
The maximum of the specific heat is usually smaller in 
restricted geometries than in the bulk, for all the bound- 
ary conditions imposed on the walls of the container, with 
the exception (the only one known by the present author) 
of Dirichlet boundary conditions . The theoreti- 

cal investigation of BEC in harmonic traps (see Ref. fl^ j 
and references therein) was motivated recently by its re- 
alization in ultracold trapped atomic gases fl3|| . In this 
situation, the specific heat of an infinite system presents a 
discontinuity at the condensation temperature. In finite 
systems, the discontinuity is again rounded off, as shown 
by analytical and numerical calculations, for example in 
Refs. [p^l7t. As the number of particles is decreased 
the condensation temperature decreases Jl4[-p7[. Fur- 
thermore, the multiple-step BEC was introduced in Refs. 

for the cases of very anisotropic boxes or confining 
potentials. In this case a finite Bose gas is condensing 
gradually to the ground state, exhibiting in between 2D 
and/or ID macroscopic populations. 

In a very anisotropic Bose system, as particle density 
decreases, the multiple-step BEC (MSBEC) temperature 
becomes lower than the temperature at which some of the 
degrees of freedom of our system freeze out. During both 
processes (MSBEC and freezing) the 3D particle distri- 
bution transforms gradually into a lower dimensional dis- 
tribution. On the other hand, the two processes change 
into each other at the variation of the particle density 
or of the dimensions of the system. Moreover, the re- 
duction of the dimensionality of the particle distribution 
due to the freezing out of some of the degrees of freedom 
can happen also for fermions at low densities. The analo- 
gies and differences between the two processes mentioned 
above, justify (arguably, of course) the use of the simpler 
expression of Bose-like condensation for the freesing out 
of degrees of freedom, in the limit of low particle density. 

The temperature at which BLC occurs (as in the case 
of BEC in finite systems, this temperature cannot be 



uniquely defined) depends on the energy spectrum and 
has a finite positive value. This type of condensation 
is identical for both bosons and fermions (see Fig. 1). 
To show this, let us consider a closed system of massive 
bosons and fermions described by a single-particle Hamil- 
tonian of the form H = H c + H^, with the eigenvalues 
e — e c + (y Ll as explained in the introduction. The mean 
occupation numbers of single particle energy levels e, are 
(rie ) = (exp (q + e/k^T) ± l) -1 , where (— ) is the su- 
perscript for bosons, and (+) for fermions, a = —[i/k^T, 
and [i is the chemical potential. We introduce the func- 
tions 



^ = E 



knT 



(1) 



da 



in a similar way as Pathria and Greenspoon did for 
bosons in j6|. Then, for example, the number of par- 
ticles, the internal energy, and the heat capacity can 
be written as = Z { ±] , U (± 1 = k B Tz{ 



(±) 



and 



C ( v ±} = k B {G^ ] - g[ ±] /G^ ±} ), respectively (in all this 
paper we shall consider spinless particles). To avoid 
divergent terms that occur in the functions introduced 
when T approaches zero, in the case when the ground 
state energy eo is positive, we redefine a as a — eo/kT and 
e as e — eo • Making these replacements we do not change 
the thermodynamics of the canonical ensemble [Q. If 
the density of the energy levels of the (quasi)continuous 
spectrum, as a function of energy, is cr(e c ), then we can 
write 



oo 

z o ±] = E 

i=0 



a(e) 



exp (a + l3ei +/3e)±l 



de , 



(3) 



where (3 = 1/fesT. If in the temperature range of interest 
for the study of BLC (ei/k^T « 1) a > 1, then we can 
write Zq in terms of two functions, corresponding to the 
continuous and to the discrete spectra, respectively: 

Z {±) - e- a Z {±) Z {±) 

where Z^ ±] = f™ a{e)e~^ de and zf> = YZo^- 
Within this approximation is no difference between 
bosons and fermions and, according to Eq. 
Zn ■ Using the relation 



r (±) 



da n 



= f3 n 



d n Z t 



(±)' 



d(3 n 



(4) 



that holds for bosons ||, as well as for fermions, we can 

write the specific heat Cy^ = Cy /N^\ in units of fee, 
as 



2 



Kb 



d(3 2 log (A a) 



(•5) 



/3 2 ^log^(/3)+/? 2 ^logZ d (/3) 



(where we have droped the superscript (±) as insignifi- 
cant in this case) . Acording to Eq. (m , the specific heat 
is nothing else then the sum of the heat capacities of two 
systems, each of them containing a single particle under 
canonical conditions, and it is described by the Hamilto- 
nian H c and H^, respectively. 
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FIG. 1. Specific heat in units of fca vs. temperature for 
ideal Bose and Fermi gases trapped inside cuboidal boxes with 
Dirichlet ((a) and (b)), Neumann ((c) and (d)), and periodic 
((e) and (f)) boundary conditions on the walls. The Figs, 
(a), (c), and (e) correspond to h — > oo, h = 10 -9 m, and 
Z3 = 10 -10 m, while the Figs, (b), (d), and (f) correspond to 
h —* 00 and h ~ h ~ 10~ 9 m. The particle density is 10 25 
m -3 in each case. In each situation the results for bosons 
(solid line) and fermions (dashed line) are both plotted, but 
they can not be distinguished. In the insets of (a), (c), and 
(e) we show low temperature details of the larger graphs (the 
axes are the same). 



Explicit expressions for Zn and C>4 can be obtained 
if we assume that the density of states of the continuous 
spectrum has the form <r(e c ) = Ce s c (C and s are con- 
stants, such that C > and s > —1), as it happens in 



most of the cases Using the Eqs. (||), and 
we can write: 



7 (±) _ 



and 

r (±) 



5sgrEi=o3ir( a + l+n-i) 



^E; = „Qr( s + i + n-j) 



(6) 



where rtj is the degeneracy of the level with energy and 
Ci = n\/j\(n — The functions (a) are the / th or- 
der polylogarithmic functions (see for example Ref. 
and the references therein for more details) of argument 
e~ a (bosons) or — e~ a (fermions). In the case of ideal par- 
ticles inside a rectangular box of dimensions l x ^> l y ,l z , 
we can write e c = % ^/2m and e^ jj = Ti (k 2 i +k 2 j)/2m, 
where fc x , k y , and k z are the wave vectors along the x, 
y, and z axes, respectively. The mass of one particle 
is m and the discrete values of k y i and k Z j depend on 
the boundary conditions. In this case s = —1/2. If 
l^,l y ^> l z , then s = and e c = h 2 (k^ + ky)/2m, while 
€i = H 2 k% i /2m. Let us now concentrate on the BLC of 
particles inside such rectangular boxes. In Fig. 1 we can 
see the results of the exact numerical calculation of cy 
(using the formulae from Eq. (||) for Z n ^ and G'£ ± ' ) ) as 
a function of temperature, for two different kinds of ge- 
ometries and for Dirichlet (Fig. 1 (a), (b)), Neumann 
(Fig. 1 (c), (d)), and periodic (Fig. 1 (e), (f)) bound- 
ary conditions. In geometry I (see Fig. 1 (a), (c), and 
(e)) l 2 = m, ^ 3 = 10~ 10 m, and h > h, while in 

geometry II (see Fig. 1 (b), (d), and (f)) l 2 = h = 1CT 9 
m and l\ 3> 1%. To make concrete calculations we choose 
A 2 = 2nh 2 /mkBT — 10 _18 T _1 which corresponds to a 
mass of about 3 atomic mass units for all the particles 
in the systems investigated. In the figure, the results for 
bosons and fermions are indistinguishable, as expected 
for low particle densities. The choice of the dimensions 
in geometry I allow us to observe the BLC from 3D to 2D 
and, at lower temperature, from 2D to ID. We observe 
the formation of a maximum (at, let us say, temperature 
Tmax) m each of these two cases and for all boundary con- 
ditions. The height of this maximum and, in general, the 
shape of the function cy(T) around T max depend on the 
spectrum of H<\- For example, for Neumann boundary 
conditions, we observe the formation of a minimum at a 
temperature a bit higher than T max . In geometry II we 
observe the BLC from 3D to ID. In this case, the max- 
ima are more pronounced and the minima observed in 
geometry I for Neumann boundary conditions disappear. 

In Fig. 2 we plot T max /T c and c Vmax /k B vs. l 3 /l, for 
Dirichlet, Neumann, and periodic boundary conditions, 
in the cases when Z 1 ,Z 2 3> h (Fig. 2 (a), (c)) and l\ ^> 
h = h (Fig. 2 (b), (d)). T c is the bulk BEC tempera- 
ture, given by the equation p(2-Kh 2 /tokbT c ) 3 / 2 = £(3/2), 
p is the particle density, £ is the Riemann zeta func- 
tion, and I = p- 1 / 3 is the mean interparticle distance. 
Since T max converges to a finite value and T c — ► when 
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p — » 0, the ratio T max /T c diverges in this limit. As p 
increases, T c increases and BLC is gradually replaced by 
BEC. As a consequence, lim /3 T max /T c = 1. Fig. 
2 (c) can make the connection between these numerical 
calculations and the analytical approximations reported 
in Ref. ||. We observe that cv max is higher in Fig. 2 
(d) then in Fig. 2 (c), at the same value of h/l, for any 
boundary conditions. Nevertheless, the maximum value 
of cvmax! which is about 2.02/cb, is obtained for periodic 
boundary conditions in the limit p — > 0, while at higher 
densities this decreases under its bulk value, as expected 
from previous calculations B. 




FIG. 2. The temperature of cvmax of Bose gases, scaled by 
the bulk critical temperature T c (see the text), as a function 
of I3/I (where I — p~ 1//3 , and p is the density), is shown for (a) 
the membrane geometry (Zi, h S> h) and (b) the wire geome- 
try (h S> h = ^3). The value of cv m ax, in units of vs. h/l, 
is plotted for (c) membrane geometry and (d) wire geometry 
(the same as in (a) and (b)). Solid, dashed and dotted lines 
are used for Dirichlet, periodic and Neumann boundary con- 
ditions. The thick horizontal lines in (c) and (d) correspond 
to the 3D bulk value of cv at the BEC temperature. In the 
numerical calculations we varied p, keeping I3 — 10 -9 m. 



The study the BLC of ideal particles in harmonic 
traps is easier since in this situation Z& has a very 
simple analytical expression. If we denote the char- 
acteristic frequencies of the harmonic trap by u x , uj y , 
and uj z , with u x <C w y ,ct> z , then Z c — k^,T '/huj x and 
Z d = [(1 - exp(^ x /A:BT))(l - exp^y/fceT))]- 1 . In 
this case dcy/dT > for any temperature, so BLC 
is not accompanied by the formation of a maximum. 
The dimensionality of the system (say, nD) is reflected 
in the value of cy, which is nfce, and the fraction of 
the particle number in the ID gs, has the expression 
N 1D /N = (1 - e -^/ feBT ))(l - e -(WkBT))_ 



III. DIVERGENT BEHAVIOR OF CV IN 
FERMIONIC SYSTEMS 

In this section we shall concentrate on Fermi systems 
close to T — K. We consider again that the Hamilto- 
nian of the system can be approximated by single-particle 
operators of the form H = H c + H^, as explained in the 
introduction. At the increase of the particle density or 
of the density of the eigenvalues of the operator 
we would expect to approach the limit in which both, 
H c and have continuous spectra (3D bulk limit). In 
such a limit we should recover the results from Ref. 0, 
namely cy a T at low temperatures. As it will be shown 
next, this is not the case in general. The continuous limit 
is not attained in a smooth way. Instead, in some situ- 
ations, the specific heat would become divergent at zero 
temperature, for certain values of the Fermi energy. 

At temperatures close to K the chemical potential of 
a Fermi system approaches the Fermi energy ep ■ For a <C 
— 1, the polylogarithmic functions of negative argument 
can be written in the form |Q : 



r(. + i)L 1 + °UJJ • (?) 

The cases for n =0 and 1 are included in (0), but can be 
refined further to write g ( n ] {a) = \a\ n [1 + (e Q )]. In 
the other extreme case, when a ^S> 1, all the polyloga- 
rithmic functions have a behavior of the form (a) = 
e~ a [1 + O (e~ a )]. Using these asymptotic expressions 
we can return to the study of the specific heat close to 
zero temperature, for a density of energy levels of H c sim- 
ilar to the one introduced in the previous section, namely 
cr(e c ) = Ce s c . The ground state of is nondegenerate 
since we discuss a finite system. We shall use the notation 

Since we know that p — > ep as T — > 0, let us now 
calculate limr^odaol — |a|) when ep = £i, i > (in all 
the other cases will turn out that the limit is zero) . Using 

N = Zy + \a), Eqs. (||), and the definition of the Fermi 
energy, we write two different expressions for the total 
number of particles in the system: 



C 



N = {s + jyi {H s+1 + ■ • • + (H " ^-i) s+1 ) 



i + o[ — 

a* 



(8) 



+n I C7^±il^ ) 1 (a + /3e J ) 



C 



j=i+i 



= (s + 1)^+1 {|ao|s+1 + • • • + (w - ^) s+1 ) ■ w 

If we denote £ = a+(3ei, then from (||) and (||), neglecting 
the exponentials and assuming that limT-^o (£/!«() I) = 
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UmT^o^/M) = Oj we obtain, in the case cto,a <C 
an equation for £: 



r(s + i 



"X-s 



(10) 



where % s = 1 + . . . + rii_i(l — Xi-i) s and Xj = tj/tF- We 
now notice that we have three distinct situations: (a) 
s > 0, in which case £ — > as T — > 0, (6) s=0, and £ 
converges to a finite positive value, and (c) s e (—1,0), 
when £ — > oo as T — * 0. 

Let us now analyze the asymptotic behavior of £ in the 
case, (c). For ^> 1 we can write 



£ n,T(s + 1) 



(11) 



so£ = (-s)log|ao|-log^-log(x a /nir(s + l)). There- 
fore, at a < -1, £ « |s| log jo;o| - log [log |a |] + . . .. 
We can see now that the assumption lim:r-»o(£/|c*o|) = 
limT^o(£/|o ! |) = was justified. Also, following the same 
kind of reasoning, one can prove that when ep 7^ £j, for 
any i, then limx->o(lQol — |a|) = for any s. 

Using the Eqs. (JlOj) and (|llj) we can calculate the 



specific heat close to K. For that we have to evaluate 



the functions G { 2 +) , g[ +) , 



r {+) 



and . We analyze 
again the case when ep = £j, i > 0. After some algebra 
and dropping out the factors that become exponentially 
small in the limit T — > 0, we can write: 



G 



(+) _ 



|s+2 



3s+l 



1 o I 4 



(12) 



TV 1 n\S*+2(£) 



(+) 



+r( s + l)yf 

C\a\ s+1 



2 9 l s + \0 



2T(s 

St 



2)y 



T, 



s+l 



r( s + 2) 



|o|" 

1 + 

51:1(0 



1 



(13) 



a 



s+l 



+ r(s + l)y 



G 



(+) _ 



l«l 
G|a| 



s+l 



(14) 



+nir(s + 1) 



C\a 



s+l 



(n + 




(15) 



where j/j = /3ej/\a\ and T s = ^l=i n fe(l _ ^fe)* 1;E fe- To 
see the asymptotic behavior, we calculate cy separately 



for the cases (a), (6), and (c). Using Eqs. ( |10] , |12 - 15) and 
working consistently in the orders of \a\, we obtain the 
following asymptotic results: 



Case (a) 



cy 



(s + l)\a\ 



k B Xs+1 +0(\a\-(°+V) 



(16) 



n^(s + l)g { s +) (0 



2s 



nfr 3 ( s + i) 3 ( +) (0 , Q ( I 



Xs l«l 3s 



a 



where m = min {s + 1, 4s, 2}. 
Cast (b) 



cy 



fc B HCxi + lodol- 1 )) 



(17) 



Xo + n».9o +) (0 Xo + n z g { Q +) (£) 



23^ +) (0 



^,9i +)2 (g) 
Xo + ^5o +) (0 



+ yX0 + O(e Q ) 



• Case ('cj 

cy 
fe B 



s + l 



Xs+i+0(|a|-^+D) ^ 



(18) 



So, for ep = e%, * > 0, from Eqs. (|16|-|18[) we distinguish 
the following situations: 

(al) s > 1/2, then cv/^b oc (ef / ksT) 1 ^ 28 , so cy — > as 
T — > (note that if s > 1 some of the orders of a 
interchange, but the function cy converges fast to 
zero as T approaches OK); 



(a2) s = 1/2, then liniT^o(cy/£:B) 
2V2)(3V8)C 2 (l/2)nf/x 3 /2Xi/ 2 ; 



(3 



(a3) s e (0,1/2), then c v /£;b oc {t F / k B T) x ~ 2s , so c v 
00 as T — » 0; 

(6) s = 0, then cy/k^, oc kBT/eF, so cy — * as T — * 0; 

(c) s € (-1,0), then cy/fce oc (e K /fc B T)/ log (e P /fc B T), 
so cy — ► 00 as T — > 0. 
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Therefore, in the cases (a3) and (c), cy presents a diver- 
gent behavior at T — K, while in case (a2) approaches 
a finite limit. These situations seem to be in contradic- 
tion with the third law of thermodynamics. To clarify 
this we mention that the divergency appears just if the 
spectrum of H c is continuous. In any finite system this is 
not the case, so at low enough temperatures cy decreases 
towards zero. 

Without getting into details we state that when ep ^ 
e,, Vi > 0, similar calculations lead us to the results 
limj^o(l a ol — M) = an d lim^^o c v = for any s. 
Moreover, in the low temperature limit we reobtain the 
known result |7| cy oc T. On the other hand, the con- 
tinuity of fx as a function of £p implies the continuity of 
a and cy as functions of €f, for any T > K. In other 
words, the divergent behavior in the cases (a3) and (c) 
can be approached asymptotically for any T > K, as 
£f — > e« (for any i), by the functions cy(T). This leads to 
the formation of a maximum at finite temperature, with 
the properties: cy max — > oo and T max — > 0, as ep — * e*, 
for any i > 1. 
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FIG. 3. The specific heat (in units of /cb) of a Fermi gas 
trapped inside a cuboidal box (Zi — » oo, h = h = 10 -9 m) 
with Neumann boundary conditions on the walls. The four 
curves correspond to the following densities: 1.5 x 10 26 m -3 
(dotted line), 9.2 x 10 26 m -3 (dashed line), 9.6 x 10 26 m -3 
(solid line), and 1 x 10 27 m -3 (thick solid line). This last case 
correspond to ef = ei. 



Let us now make the connections with familiar sys- 
tems, namely with the ones discussed in section |n|. In 
the case of a cuboidal box with dimensions l x 3> l y ,l z , 
s = —1/2, so we are in the case (c). In Fig. 3 we plot 
the exact numerical calculation of such a fermionic sys- 
tem, with dimensions h — > oo, li — h — 10~ 9 m. The 
mass of the particles is chosen as in Section [H], such that 
A 2 = 10 -18 T _1 . We observe the formation of the max- 
imum as the Fermi energy approaches the first excited 
energy level of H^, and the divergent behavior at ep = €\. 
If the fermions are inside a cuboidal box with dimensions 
l x , ly S> 1 2. or a harmonic potential with the characteristic 



frequencies u> x <C u> y ,uj z , then s = and we are in the 
case (b), therefore we do not observe the formation of a 
similar maximum. This was checked by exact numerical 
calculations and was found to be correct. 



IV. CONCLUSIONS 

In Section [H] of this paper it is presented in general 
the phenomenon of Bose-like condensation in the case of 
massive bosons and fermions. This denomination was in- 
troduced in H where, according to my knowledge, it was 
reported for the first time a crossover between different 
dimensionalities of the phonon gas distribution in ultra- 
thin dielectric membranes. This phenomenon appears to 
be identical for both types of massive particles and re- 
sembles to the multiple-step Bose-Einstein condensation 
HQ- Nevertheless, the two phenomena are different in 
nature. The results are exemplified for the familiar cases 
of ideal particles trapped inside cuboidal boxes and har- 
monic potentials. 



The analysis made in Section III, lead us to the ob- 



servation of interesting divergences of the specific heat of 
a Fermi system at zero temperature. The phenomenon 
is described in general, for a single-particle hamilto- 
nian of the form H = H c + H^, with H c having a 
(quasi)continuous spectrum e c € [0, 00) with the energy 
levels density cr(e c ) = Ce s c (s > —1) and having the 
discrete eigenvalues ej, i = 0, 1, . . .. It was found that 
c v (T) -> 00 as T -> for any s G (-1,0) U (0,1/2) 
if ep = ej, for any i > 1. This divergent behav- 
ior is approached asymptotically for any T > 0, as 
£f Vi > 1, leading in this way to the formation 

of very high maxima (in the limit, infinitely high) of the 
fermionic specific heat close to zero temperature. This is 
an unexpected new phenomenon, since it seems to con- 
tradict the third law of thermodynamics. Anyway, this 
does not happen since in any finite system the energy 
spectrum is discrete and at low enough temperature the 
specific heat decreases towards zero. Nevertheless, this 
phenomenon might have interesting consequences on the 
entropy of the system in the vicinity of zero temperature. 
On the other hand it should be investigated if systems 
obeying fractional-statistics ]22] o r interacting Bose sys- 
tems (see for example Ref. p3| , |24| ] and references therein 
for similarities between these two types of systems) ex- 
hibit similar behavior. 
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